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The analysis of the size distribution of droplets condensing on a substrate (breath figures) is a 
test ground for scaling theories. Here, we show that a faithful description of these distributions 
must explicitly deal with the growth mechanisms of the droplets. This finding establishes a gate- 
way connecting nucleation and growth of the smallest droplets on surfaces to gross features of the 
evolution of the droplet size distribution. 
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Classical questions regarding breath figures involve the 
influence of material defects and impurities on the droplet 
patterns [1-5]. Presently, they are used as self-assembling 
templates in micro-fabrication [4, 6-10], as highly effi- 
cient means for heat exchange in cooling systems [5, 11- 
13], and they are promising candidates for water recov- 
ery in (semi-)arid regions [14, 15]. For these applica- 
tions a detailed knowledge of the droplet size distribution 
and the average droplet growth speed is vital. Here, we 
demonstrate that the state of the art scaling theory [16- 
20] fails to describe data from simulations [22] and lab- 
oratory experiments [21], Fig. 1. A faithful description 
must therefore explicitly address the microscopic growth 
mechanisms of droplets. 

Classical scaling [16, 17] asserts that on clean surfaces 
the coagulation of droplets organizes the systems into 
a state where the number of droplets, n(s,t), per unit 
droplet volume and surface area takes a universal scaling 
form, 

n(s,t) = S - e f(^) , with S = S(t). (la) 

Here s denotes the droplet volume, 9 is a scaling expo- 
nent, f(x) is a dimcnsionlcss function, and S(t) is the vol- 
ume of the largest droplets encountered at time t, i.e., the 
average volume of droplets in the bump of the distribu- 
tions shown in Fig. 1. 

Since n(s, t) has a dimension of length to the power —5 
the exponent 9 must be set to a value of 9 = 5/3 [17-19]. 
The time evolution of S(t) is found by observing that 
the total volume of all droplets grows linearly in time 
when a constant volume flux impinges onto the surface. 
In agreement with experimental and numerical observa- 
tion [17-19] this entails S(t) <~ t 3 . Moreover, a lower 
cutoff to the scaling at a scale sq/S has been accounted 
for by a polydispersity exponent < r < 2 [23]. For our 
numerical scheme, where the mass flux onto the surface 
is implemented as sustained addition of droplets of size 
so to random positions of the surface and where overlap- 
ping droplets are subsequently merged [21], it was pre- 



dicted [20] to be 

x < 1 f(x) ~ x e - T , t = 19/12. (lb) 

The scaling, Eq. (la), provides an excellent data col- 
lapse of the bump and the dip of the numerical, Fig. 1(a), 
and the experimental data, Fig. 1(b). Beyond the dip one 
can discern a self-similar scaling regime, Eq. (lb), in the 
numerical data, and in only those experimental data with 
the vastest range of droplet sizes. On the other hand, 
in either case — and particularly pronounced in the ex- 
perimental data — noticeable deviations, tails, from the 
scaling prediction arise for small values of s/S. 

In the following, we show that these deviations result 
from features of droplet growth at the small length scale, 
so- Similar to the approaches in the theory of critical 
phenomena [24, 25] or of the effect of rough boundaries 
in turbulent flows [26, 27], scaling will be recovered by 
asymptotic analysis [28] , which allows us to explicitly ac- 
count for different growth mechanisms of small droplets. 
Universal and non-universal features of the asymptotic 
droplet density distribution will be disentangled by dis- 
cussing the consequences of the different growth mecha- 
nisms for the small droplets in the numerical and exper- 
imental setting, respectively. 

Relation to fractal packings. — To explore the role of 
the lower cutoff of scaling we consider the droplet ar- 
rangement in breath figures as an example of a fractal 
packing of disks (see [29-31] for recent applications in 
other fields), and adopt scaling arguments developed to 
characterize (disordered) fractal structures to the prob- 
lem at hand: We assert that in the scale-separation limit, 
so <C S, the free surface area, i.e., the area not covered 
by droplets, approaches a fractal with a fractal dimen- 
sion df < 2. Self-similarity with fractal dimension dj 
amounts then to the statement that in an area of size 
5 2 / 3 a number 
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FIG. 1. (color online) Scaling plots of the droplet number density n(s,t) for (a) numerical and (b) experimental data, for 
(a) 8 x 10 6 (rightmost tails), 12 x 10 6 , 16 x 10 6 , 20 x 10 6 , and 24 x 10 6 (leftmost tails) droplets added to a domain of size 
1600 x 1600 [21], and (b) water droplets on a polyethylene film [22] where eventually the droplet diameters cover the range 
from a few microns to a few millimeter. The respective positions of the tail, dip and bump of the distributions are indicated, 
and the slope, 9 — t = 1/12, of the scaling prediction, Eq. (lb), is shown by solid lines. The insets show snapshots of the breath 
figures at an early [blue, (lower) right] and a late [red, left] time, respectively. [MOVIES] of their time evolution and full details 
of data assimilation and evaluation, as well as plots of the raw data are given in [SUPPLEMENTARY MATERIAL] . 



of regions of size s 2 / 3 are required to cover the comple- 
ment of the surface area covered by all droplets larger 
than s*. The fraction of this area in the considered do- 
main of size S 2 ' 3 amounts to 



p(s*,S)=N(s*,S) 
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Following [29] we denote the surface area not covered 
by droplets as porosity, p(t). It is obtained by evaluating 
Eq. (3) for the size Sq characterizing the small scale cutoff 
of the fractal, p(t) = p(sq, S(t)). 

By its definition the porosity is related the area Ad 
covered by droplets in a region of area A s via p(t) — 
1 — Ad/A s : when the surface area in between droplets 
approaches a fractal of zero measure one obtains 



a{s) n{s, t) ds = — - = 1 — p(t) > 1 . 

Af, 



where a(s) denotes the area covered by droplets of size s. 
Using Eq. (1) with = 5/3, Eq. (lb), a(s) ~ s 2 / 3 , and 
introducing x = s/S one obtains 
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and comparing Eqs. (4) and (3) yields 



(2-dy)/3. 



(4) 



(5) 



Hence, the nontrivial scaling, Eq. (lb), of f{x) for small 
x reflects the fractality of the arrangements of droplets 



in breath figures with a large scale separation s <C S, 
that can faithfully be regarded as a fractal. 

This provides an independent, more accurate, means 
to test the polydispersity exponent: For t = 19/12 and 
S ~ t 3 , Eq. (4) implies p(t) ~ i -1 / 4 . Remarkably, none of 
our data follow this prediction [Figs. 1(b) and 5(b) in the 
supplementary material]: Rather than 1/4 we find 0.30 
for the numerical, and 0.16 for the experimental data. 

Hence, the different microphysics of droplet growth 
and merging leads to (slightly) different fractal dimen- 
sions and a different small-scale cutoff of scaling. To 
disentangle the intermediate self-similar scaling regime 
from the large scale (arising from the first generation of 
droplets, cf. [32]) and the small scale physics we intro- 
duce cutoff functions f(s/S) and g(s/ s o) f° r large and 
small droplets, respectively: f(x) = f(x)/x e ~ T takes a 
constant value /o for x <C 1, and it accounts for the dip 
and the bump in n(s,t) for s ~ S. Similarly, g(s/so) ac- 
counts for the tails of n(s,t). As shown in the insets of 
Fig. 2 it approaches constant values for s > so, and it 
takes a scaling form for all times. To arrive at a complete 
description of the droplet size distribution we further dis- 
cuss now this lower cutoff. 

Incorporating the lower cutoff. — We start by writing 
the total volume of droplets of size s per unit volume ds 
and unit surface area in the form 

s n(s, t) = S- 2 / 3 (|) ~ d,/3 f(s/S) g(s/s ) (6) 

This expression states that in the scaling regime so "C 
s«S the overall volume of droplets of size s on an area 
of size S* 2 / 3 is proportional to the number of droplets, 
Eq. (2), of the considered size. 

In equation (6) the fractal dimension df and the func- 



3 



(a) 



10 u 



to 



10" 



10" 




10 10 10 10 

s [mm 3 ] 



10" 



10" 4 10" 

s/S 



10° 



FIG. 2. (color online) Master plot of the universal cutoff 
function, f(s/S), which comprises all numerical and exper- 
imental data shown in Fig. 1 (using the same symbols), as 
well as additional numerical data shown as black crosses (ev- 
ery 5 x 10 4 frames of 44 runs where 4.23 x 10 7 droplets are 
added eventually). The insets show g(s) for (b) numerical and 
(c) experimental data at all times. Their shapes (solid lines) 
reflect the respective small-scale droplet growth mechanisms, 
Eqs. (10) and (11). The vertical line in (c) marks so- 



tion g(s/so) are not universal. We henceforth adopt the 
values for df determined by fitting the porosity, and we 
follow the evolution of droplets smaller than over a 
small time interval from t to t + dt in order to relate the 
form of g(s/ Sq) to different nucleation and growth mech- 
anisms of small droplets. Let the considered area on the 
substrate and the time interval dt be chosen such that 
this growth is not influenced by the merging of droplets 
larger than s*. When s* is so small that f(s*/S) takes 
the constant value /o, Eq. (6) entails 



n 



In order to determine <?(s*/so) we observe that the vol- 
ume density of the droplets smaller than amounts to 
the cumulative distribution 



V(a m ) 



s n(s, t) ds 
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-df/3 



j(s/s )ds. (7) 



For values of in the scaling regime the increase of vol- 
ume is accounted for by increasing the integral domain. 
Therefore, an infinitesimal increase of s* to s* + s dt in 
the time interval dt amounts to an increase of the volume 
density of droplets, 

dV _ V(s*+sdt)-V(s*) 
dt dt 



fo 
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On the other hand, this change must be due to the vol- 
ume flux $ onto the fraction of area covered by droplets, 



dV 
~dt 



K*.,s)$~(£) 



(2-d / )/3 
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Equating the expressions for dV/dt, Eqs. (8), and drop- 
ping the subscript *, one hence obtains 



g(s/s ) 



$ s 2/3 

To * 



(9) 



This expression provides the desired connection of the 
average speed s of the growth of droplets of size s to 
the form of the small-scale cutoff g(s/so) of n(s, t): it ex- 
plains how different microscopic droplet growth laws give 
rise to different non-universal cutoff functions g(s/so), 
and how the universal scaling is recovered for s 3> sq. 
After all, the volume growth of large droplets in breath 
figures is always proportional to the area exposed to the 
surface flux [16, 17, 20], s ~ $ s 2/3 . Equation (9) allows 
us to disentangle universal and non-universal contribu- 
tions to n(s,t). This major finding of our theoretical 
treatment is now substantiated by working out the mul- 
tiscaling predictions for the data shown in Fig. 1. 

Scaling numerical data. — When a small droplet, of 
size so j is been added to the surface, it is merged with a 
droplet on the surface when the droplets overlap. As a 
consequence, a droplet of radius s 1 / 3 will capture small 

1/3 

droplets of radius s , that are added in a distance 



smaller than s 1 / 3 



1/3 



from its center. In the absence 
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of other droplets, this growth amounts to 



The term in square brackets accounts for an enhanced 
growth of small droplets s> s , which ceases rapidly for 
increasing s. In practice the decay is even faster since 
the capture regions of neighboring droplets overlap. To 
fit the simulation data, Fig. 2(b), one therefore needs a 
non-trivial prefactor 0.76 and an exponent close to 0.78 
rather than 1/3, 



g(s/s ) ~ 0.07 



0.78 



(10) 



Using Eq. (6) and 9 — t = 0.3 this provides a perfect data 
collapse of all numerical data, Fig. 2(a). 

Scaling experimental data. — In the experimental set- 
ting, the growth rate of the droplets has two contribu- 
tions. For small droplets, growth is limited by the diffu- 
sion of water molecules on the substrate towards the con- 
tact line of the droplet. As derived in [33] and observed 
in the experiments of [32], the radius of small droplets 
grows then like r ~ t 1 / 4 , such that s ~ $ s -1 / 3 . For 
larger droplets, the volume flux from the vapor phase 
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onto the droplets is again proportional to the exposed 
droplet surface, such that s ~ $ s 2 / 3 . These growth con- 
tributions combine to 



$ s 2 / 3 (l + 



g(s/s ) 



(ii) 

where So — 1.5 x 10 6 mm 3 is the crossover size scale 
and b ~ 2 x 10~ 2 is a normalization constant. Inserting 
Eq. (11) into Eq. (9) provides an excellent prediction for 
g(s/so), Fig. 2(c). Also for the experimentally measured 
droplet size distributions one hence obtains a perfect data 
collapse of the appropriately scaled droplet number den- 
sity n(s,t) for all different times, Fig. 2(a). 

Discussion. — For the numerical data, s amounts 
to the volume of the smallest droplets in the system, 
Fig. 2(b), and for the experimental data it is about one 
order of magnitude larger than the smallest observed 
droplets, Fig. 2(c). In either case g(s) saturates for 
s > 10 2 so. On the other hand the scaling behavior 
Eq. (lb) is only accessible for values of s below the dip of 
the distribution, i.e., for s < 10~ 2 S. It can hence only be 
resolved in simulations where 1CT 4 <C So/S, resulting in 
the observed scaling regime of about 1-2 decades, in the 
numerical data, Fig. 1(a), and in the experimental data 
with the largest accessible scale separation, Fig. 1(b). 

Due to the relatively small scaling range the droplet 
size distribution of breath figures can not merely be ideal- 
ized as a self-similar process with a single relevant length 
scale S(t) [16-19]. Rather one explicitly has to cope with 
the growth law of the smallest droplets in the system. Via 
its (slight) effect on the fractal dimension characterizing 
the free space in between the droplets, Eq. (5), it sets the 
value of the polydispersity exponent r, and it leads to 
massively different cutoff functions g(s/s ), Fig. 2(b,c), 
that can completely dominate the shape of the droplet 
size distribution, Fig. 1(b). 

When both the large scale and the small scale cutoffs 
are properly accounted for via Eqs. (10) and (11), a re- 
markable data collapse of all experimental and numerical 
data into a single plot is achieved, Fig. 2(a). This re- 
covery of scaling, establishes a novel gateway connecting 
features of the microscopic droplet growth on surfaces to 
gross features of the evolution of the droplet size distri- 
bution. 
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